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^v^j , We consider the static wall approximation to the dynamics of a particle bouncing on a periodically 
oscillating infinitely heavy plate while subject to a potential force. We assume the case of a potential 
^ given by a power of the particle's height and sinusoidal motions of the plate. We find that for powers 
smaller than 1 the set of escaping orbits has full Hausdorff dimension for all motions and obtain 

f*^j existence of elliptic island of period 2 for arbitrarily high energies for a full-measure set of motions. 

«^ Moreover we obtain conditions on the potential to ensure that the total (Lebesgue) measure of 
elliptic islands of period 2 is either finite or infinite. 

q 

^ . I. INTRODUCTION 

a 

This work provides some results about the description of trajectories of a parti- 
cle in a model that presents Fermi acceleration. This acceleration mechanism was 
originally proposed in 1949 by Enrico Fermi [Fe49] to provide an explanation for 
presence in the universe of high energy particles known as cosmic rays. Such par- 
ticles are believed to gain energy by consecutive interactions with irregularities in 
an otherwise stationary magnetic field. Naively one would expect that a process 
of thermalization would occur, leading to a stationary motion of the particle itself. 
However, this reasoning turns out to be too simplistic and a more refined analysis 
shows that there is a definite probability of an average gain in energy. 
In 1961, Ulam [U1G1] suggested a simple Hamiltonian system to model such statisti- 
cal acceleration behaviour. The model has been thereafter known as the Fermi- Ulam 
ping-pong model and consists of a particle moving between two infinite walls that 
are performing an oscillatory motion; the particle changes its velocity only by elastic 
collisions with the moving walls and it is not subject to any other force. The main 
questions about this problem regarded the existence of trajectories with unbounded 
energy: such orbits can either be escaping, i.e. such that the energy of the parti- 
cle goes to infinity with time, or oscillating, meaning that the limsup of energy is 
infinite while the liminf remains bounded. In 1977, KAM theory has been used to 
provide ([Pu77], [Pu95]) the answers to such questions: for sufficiently smooth mo- 
tions of the wall, all orbits are bounded, as there are invariant tori for high energies 
that prevent diffusion. It is interesting to note that the smoothness condition is not 
a mere technical issue as for less regular motions one can find (many) unbounded 
orbits; see for instance [Zh97]. 

A variation on the same theme involves a single oscillating wall and a potential 
U (x) = x a to bring the particle back to the wall. This model allows for more flex- 
ibility as we can vary the potential that acts on the particle. The case of gravity 
potential (a = 1) has been investigated [Pu77] and the study yielded the following 
result: 

Theorem 1.1 (Pustylnikov). There is an open set of wall motions <p(t) (in the 
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space of periodic analytic functions admitting an analytic continuation to a given 
strip \Qz\ < e) such that the measure of the escaping orbits is infinite. 

The case of elastic potential (a = 2) has been studied as well in [Or99], [Or02]; 
abundance of unbounded escaping orbits has been proved under some resonance 
condition between the motion of the wall and the potential. 

In a more general setting we can again use KAM theory [Do07] to prove the following 
result: 

Theorem 1.2 (Dolgopyat). If a > 1 but a ^ 2 and the motion of the wall is 
smooth enough, then the set of escaping and oscillatory motions is empty as there 
is persistence of invariant tori for high energies. 

On the other hand, orbits with unbounded energy are allowed for weaker poten- 
tials. It is conjectured that for potentials weaker than gravity (i.e. for a < 1) the 
measure of escaping motions is zero. A first step towards this result is the following 

Theorem 1.3 (Dolgopyat). If a < 1/3 and the motion of the wall is a sinusoid, 
then the set of escaping orbits has zero Lebesgue measure. 

The above results leave several open questions regarding the largeness of the 
following sets: 

• The escaping set S i.e. the set of orbits such that energy tends to infinity as 
time grows; 

• The set of orbits with bounded energy; 

• The oscillatory set i.e. the set of orbits such that limsup-E'(f) = +oo and 
liminf E(t) < oo. 

In this paper we consider a static wall approximation to the bouncing ball system. 
This approximation, described in more detail in section II, is widely used in physics 
literature. It has the advantage of being given by more simple and explicit formulae 
whereas keeping all essential geometrical features of the complete model. The first 
result of this paper is the following 

Theorem A. If a < 1 and the motion of the plate is a sinusoid, then the escaping 
set § has full Hausdorff dimension. 

As a matter of fact we still ignore what should be the measure of the escaping set 
$ in the case 1/3 < a < 1. The proof for the a < 1/3 case proceeds by defining a so- 
called critical set C of finite mass for some invariant measure, and such that almost 
every orbit passes infinitely many times through C; applying Poincare recurrence 
theorem to the first return map to C gives the desired result. The key ingredient 
of the proof is that orbits outside C are strongly hyperbolic. The set C has infinite 
measure for a > 1/3, therefore we aim to better understand dynamics on C in order 
to obtain a smaller set C of finite measure with the same recurrence property As a 
first step we study elliptic islands inside C; more exactly we establish some results 
about the abundance of elliptic islands for high energies: 

Theorem B. Let a < 1. Then in the static wall approximation, for almost all 
sinusoidal motions there are elliptic islands of period 2 for arbitrarily high energies. 
Moreover if2/3<a<l we can prove the same result for all sinusoidal motions. 
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As a byproduct of the construction involved in the proof of Theorem B we obtain 
conditions on a for which the total measure of elliptic islands of fixed period 2 is 
either infinite or finite. We remark that Theorem A extends without difficulties to 
the non-approximated bouncing balls system as we will briefly explain at the end 
of section III. We also expect Theorem B to remain valid for the exact bouncing 
ball system, but the computations needed to verify this statement are quite cum- 
bersome. Questions regarding existence and abundance of oscillatory orbits are still 
completely open. 



II. THE MODEL AND MAIN APPROXIMATION 



We consider the problem of a point mass bouncing vertically on an infinitely 
heavy horizontal plate which oscillates with period 1 in the vertical direction and 
interacts with the particle by the law of elastic reflection. The particle is moving in 
a potential U(x) = x a , where x is the vertical position and a is some positive real 
number. 

Let 4>(t) be the vertical position of the plate at time t, periodic of period 1; for 
simplicity we will consider the case: 

<j>(t)=B+£-sm(2wt). B>A (1) 

Z7T Z7T 

It is natural to associate to the system a discrete time map defined as follows. Let 
t n be the time of the n + 1st collision between the plate and the particle and v n its 
velocity (pointing upwards) immediately after the collision. As the position of the 
plate is a 1-periodic function of time t we can consider t n on S 1 = K/Z. In this way 
the phase space is a half cylinder, as the velocity immediately after a collision at 
time t has to satisfy the inequality v > (/>(t). 

An approximation common to this kind of problems is the so-called static wall 
approximation, in which one considers the plate fixed at position xq but exchanging 
momentum with the particle as if it were moving. Notice that since the displacement 
of the plate is a bounded quantity, it affects the map by terms of order at most 1/v, 
i.e. the accuracy of the approximation increases for high energies. 
Define T(v) : M + — > M + as the time taken by a ball leaving x = xo with upward 
velocity v to return on x = xq subject to the potential U(x). In our case T(v) ~ i> 7 , 
where it is easy to check that 7 = 2/ a — 1. In fact, let Eq be the energy of the 
particle after the collision: 



l r^** 1 

T = 2 I —dx 
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v(x) ' J xo V^o - x° 



J 



Performing the change of variable x = XQ^ a y we obtain: 

T = £ s _ ia \ _ El' a dy = B =-»Con8t • ( 1 + ( 



^-1 



Notice that the asymptotic expression is exact for xo = or a = 1. 

We will study the system in the static wall approximation and considering T = C-v 1 . 

In this setting the map F : (t n ,v n ) 1— ► (t n +i,V n +i) can be written as: 



(2) 
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Note that if we consider T{v) rather than v as the fundamental variable we obtain: 

T ^ C ((§) / +2 ^ t + T ^j ~ + T 1-1 / 7 C' <fi(t + T) + h.o.t. 

This map resembles the standard map, with the perturbation parameter k = 
rpi-ihc 1 depending on v. We can distinguish between the following cases: 

• weak potentials 0<a<l«-»7>l: the perturbation parameter grows as 
velocity grows; we expect diffusion to high energies; 

• gravity a = 1 <-> 7 = 1: the system is equivalent to the standard map 
(unfolded on a semi-cylinder along v); 

• strong potentials I<a<2«-»0<7<1: the perturbation parameter 
decreases as we increase energy; there is persistence of invariant tori for high 
energies and therefore we can not have diffusion to arbitrarily high energies; 

• elastic potential a = 2 <-» 7 = 0: the period T is independent of v, this is an 
exceptional case; 

• strong potentials (II) a > 2 <^ — 1<7<0: high energies correspond to 
small T, therefore we recover once more invariant tori bounding energies from 
above; 

• Fermi-Ulam ping pong a — > 00 <-> 7 = — 1: limit case of strong potentials; 
once more high energies are bounded by invariant tori. Notice that in this 
case the static wall approximation is not anymore a good approximation as 
we would neglect term of the same order as T. 

We will prove our results for the map (2) in the weak potentials regime. 



III. HAUSDORFF DIMENSION OF THE ESCAPING SET 



A. Main definitions 



We recall the definition of Hausdorff dimension of a metric space. First we need 
to define the Hausdorff s-measure of a metric space E as: 

H 3 (E) = lim sup inf \ V diam ( Ai) s \ . 

S^O A = {A Z } \£f 

6— covering of E ^ ^ 

Then we define the Hausdorff dimension of E as that critical s such that: 

dhn H E = inf{s s.t. H S (E) = 0} = sup{s s.t. H S (E) > 0}. 

It can be actually proved that if s < diuiH(E) then H S (E) = 00; moreover, Haus- 
dorff dimension is a biLipschitz invariant of metric spaces. 
Now let § be the set of escaping points, i.e. 

& =f {(^o,^o) s.t. v n — ► 00 as n — > 00}. 
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Theorem A. If a < 1 and 4>(t) is given by equation (1), then dim^f? = 2. 

Proof. The actual proof will be provided in the two following subsections; hereby we 
outline all necessary steps. As in [Do07] we define a set C such that orbits avoiding 
this set are hyperbolic; wc will then consider a curve T obtained as a graph of a 
function of t, transverse to the stable direction on hyperbolic points. The core of the 
proof lies in proposition III. 1 and lemma III. 8 which prove that dimjj {§ n T) = 1; 
using transversality to stable directions we obtain the desired result dim^t? = 2. 
As trivially dim# (S n T) < dim^T = 1, it suffices to obtain the reverse inequality; 
such lower bound will be found by bounding the dimension of a subset containing all 
trajectories such that the particle gain at least energy e on all (but a finite number) 
of collisions. The latter subset can in turn be expressed as the limit inferior of all 
preimages of the set of trajectories for which we gain energy at least e on the first 
collision. 

Let 7r be the projection from the phase space onto the time coordinate t and let 7 
be the parametrization of T that is a section for 7r, i.e.: 

7r : (t, v) 1— » t no 7(i) = t. 

We are going to consider the sequence of maps of S 1 to itself defined by the following 
composition: 

F„ = 7T O F n O 7. 

The proof of the result for this sequence of maps relies on the fact that on the hyper- 
bolic set the map F is expanding along T (and its images) with an expanding rate 
that increases with energy; since we are looking for points belonging to the escaping 
set, velocity will grow along with n. Therefore we need to bound the dimension 
of the intersection of all preimages of a fixed set by a sequence of expanding map 
which expansion rates are diverging. In such a setting we mimic the calculation of 
the Hausdorff dimension of the usual middle third Cantor set. Roughly speaking 
one can define a (n, A:)-Cantor set by dividing the unit interval in a number of equal 
pieces n, removing n — k of them and then iterating the procedure on the k leftover 
intervals. The dimension of the set obtained in this way is well-known to be: 

logfc 

Onk = ■ 

logn 

Notice that if n and k grow while their ratio remains bounded, the above formula for 
&nk approaches 1; the same conclusion holds true if we let n and k grow to infinity 
as we iterate the building procedure. Another way to define such Cantor sets is to 
obtain them as intersection of preimages of a given set under an expanding map of 
the interval; for example the classic middle third Cantor set might be obtained by 
taking the intersection of all preimages of the set [0, 1/3] U [2/3, 1] under the map 
x 1— > Zx (mod 1). If we let the expansion rate grow along with n wc obtain once 
more a set of Hausdorff dimension 1. To conclude the proof we bound in a biLipshitz 
fashion the sequence of maps {F n } between two sequences of linear expanding maps 
whose expansion rates are increasing, and then obtain the corresponding estimate 
for the Hausdorff dimension of the set of escaping orbits. 

In the next subsection we will prove our result for a model system; in the subsequent 
subsection we will show that the original system can be reduced to such model 
system. 
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B. Model system 

Our model system will be a sequence of continuous functions /„ : S 1 O, such that 
any lift /„ satisfies the following inequalities for all real x, y: 

m n \ x - V\ < \fn(x) - f n (y)\ < m n \x - y\, l<ra„<m„, m n ,ra„ / oo. 
Now define F. ' n = f n o • • • o f±. Then clearly we have: 

n n 

M n \x-y\<\F n (x)-F n {y)\<M n \x-y\, M n = \[m k , M n =\\m k . 

k=l k=l 

For each n we can decompose S 1 as a disjoint union of intervals of the form [a, b) 
such that F n is 1-1 and onto S 1 on each of such intervals. Each interval corresponds 
to a different choice of an inverse function for F n . In order to choose an inverse 
function for F n we have to choose an inverse branch for each one of the of which 
F n is the composition. In such a way we obtain a natural labeling for the inverse 
branches of F where we write: 

Vn S x = □ I h ... jn , 

31— jn 

and each jk runs on all different inverse branches of fk- Fix i9 £ (0,1] and an 



Figure 1: Choosing branches of Fk for k — 2. 

interval Jo C S 1 of length then define Jk ~ {Fk)^ 1 Jo — U Jji—jk where J 3 \—j k = 
Jk H Ij 1 ...j k ; finally we can define: 

J = lim inf Jk ■ 

k — 'oo 

Proposition III.l. //Vn G N m n < Cm n , then the Hausdorff dimension of J is 
1. 

Proof. In order to compute dim#J we will consider a smaller subset J'. We will 
define J' inductively as the limit of a decreasing sequence of sets J' n . First, as m n 
and m n are increasing and diverging, from some k on we have rnkfi > m k , d > 2. 
Without loss of generality we may assume k = 1. Define J' Q = J ; then, assuming 
we already defined J' n _\ and recalling that J„ = |J Jj 1 ...j n we define the decreasing 
sequence: 

Jn = I I J 31— in i 
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where the union is taken only over intervals which fully intersect J' n _\, i.e. all 
intervals that are not truncated by the intersection. In other words: each set J„ 
is the union of many intervals; some of these intervals are completely contained in 
J' n _l, and such intervals will be called good intervals and belong to what we will 
call good branches. All other intervals, i.e. such that they are not fully contained 
in J' n _i belong to bad branches. Now we define: 

J' = lim inf J' k = (| J' k . 

k— »oo 

feGN 

Clearly V k J' k C Jk ■ We only need to show that these sets are non-empty. Again 
we proceed by induction; Jq is trivially non-empty. Notice that since F~ 1 (Jo) = 
Fn-iifn 1 (Jo)) , all good intervals of J n are mapped by F n _\ to intervals of f~ x 
that are completely contained in Jo- The set Z" 1 Jo is the union of several disjoint 
intervals and since to„i? > 2, there exists at least one Ii n C Jo, therefore at least 
one of the Ji intervals will fully intersect Jo; moreover, if we let k n be the number 
of such intervals; we have the following estimate for k n : define k n: k n £ N according 
to this relation: 

k n + 1 < m n ti <k n + 2 k n + i< m n -d <k n + 2 (3) 

then we get k n < k n < k n + 2. In particular k n is always positive. 
Now to obtain J' n it suffices to take, for any such interval, the preimage relative to 
each good branch of F~_^, as J' n _i is non-empty by inductive hypothesis we have 
that J' n is non-empty as well. Moreover the number K n of intervals present in J' n 
is easily obtained by the previous inductive procedure: 

n 
1=1 



Y[ki=K n < K n < K n =Y[k t . 

i=i i=i 

With a slight abuse of notation we will write from now on J' = I I . . J,' „■ where 
now each ji is running over just the good branches. The sets J' n form a decreasing 
sequence of nonempty compact sets, therefore their intersection J' is nonempty. 
We will now show that J' C J has Hausdorff dimension s = 1. First of all it is 
obvious that s < 1 being J C S 1 , therefore it suffices to show that for all s < 1 we 
have H S (J) > 0. 

Definition III. 2. We define the running Hausdorff dimension of {ji • • ■ j n } the 
real number Sj- L ...j n such that 

The running Hausdorff dimension depends on {m r . m n } and $ according to this 
estimate: 

Lemma III. 3. Suppose that the following holds for all n € N: 

m„ < Cm n . 



Then we obtain the following lower bound for the running Hausdorff dimension: 

(4) 



log i9 + n(logf -logC) 

1 H == = 1-£ B /1 asruoo 

logM„ - log 



Proof. Using (3) we get 

k n < i)m n < k n + 2 < 3/c„ fc„ < i?m„ < k n + 2 < 3/c„. 

=K„ < tf"M„ < TK n K n < d n M n < TK n . 
Now since we know that t)/M„ < \J'j 1 ...j n |j we obtain: 

(t?/M„) Sji - 3 " < tf„ (tf/M„) s ""'" < ^„|Jj 1 ...,J s - = 1, 
and taking logarithms we establish the following inequality: 
log£„ - s h ... jn (logM„ - log i?) < 0. 
Therefore we obtain the bound: 

logM„ + ^log| 



logM n -logtf 

Now by hypothesis we know that M n < C n M n , therefore log M n < n log C+logM n 
that in turn implies: 

> logM» + n(logf -logC) 
Sjl > logM„-logtf 

that implies estimate (4) provided that we show that e„ is going to 0; but 

n 

log M n 

as log M n /n is the average of the diverging sequence logm„. □ 

To obtain a lower bound on the Hausdorff dimension we are going to use the 
following proposition (see e.g. [Fa85]): 

Proposition III. 4. Suppose there exists a probability measure fx on a metric space 
X on the a-algebra of Borel sets such that for all sufficiently small balls B we have: 

fi(B) < Cdiam(B) s , (5) 

then dimnX > s 

We will now prove the following 

Proposition III. 5. There exists a probability measure [i on J' satisfying (5) for 
all s < 1 . 
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Proof. First we build a probability measure /j on J', and then check that [i satisfies 
(5) for every s < 1. For each n and choice of j\ ■ ■ ■ j„, fix a point Xj 1 ...j n <G Jj 1 ...j n - 
Then define the following sequence of positive functionals acting on 6(§ 1 ,1R): 



•Jn, 



Now we argue that this sequence of functionals has a weak limit for n — > oo. 
In fact any continuous function 95 on S 1 is also uniformly continuous; therefore 
Ve 36 such that \x — y\ < 5 implies \tp(x) — <p(y)\ < £■ Now take n such that 
< 6. Then for each m > n: 



" lax j \J'n-j 



|*n(v) " *m(¥>)| 



1 



E E 



Jl, 



1 

Krr 



< 



< 



E 

'i,— J. 

E 

E 

ji,— .jr. 



hf^'' • 1 E ^- 

\ Jn+l,-,Jm 



< 



1 

i^ 

1 



IW^Jl-J,,) -^Ojl-jvjl 



< 



r. 



In the above inequalities we used the fact that by definition yV K n /K m 

e J' n ...j n and as | J^...^ 



and that J'j 1 ...j m C J'j 1 ...j n which implies Xj 1 ... 0n ,^ jl ... :j 



1 

( •A/ 1 -j„ ana as I^.^J < 

<5, the inequality follows. The sequence $ n weakly converges to a positive functional 
4> i.e. to a Borel measure via the Riesz representation theorem. Moreover fx is a 
probability measure (it suffices to compute the limit against the function (p = 1). 
At this point for any Borel set E and n S N we can write: 

ii{E) = lim V — XEfe...j ra ) 

Jl,— ,Jrn 

= E Um E 1T XE ( X h-3nj n +l-jm) 

Jl,— ,Jn J»+l,— ,jm 



< 



Jl,---,Jn J»+l, 



Jl,--- ,J« 



1 

if n ' 



By definition of running Hausdorff dimension this implies the following estimate for 
ME): 



ME) < £ 



I /' ■ 

rji-j„ 



(6) 



Jl,— Jn 
J il-3„ ni5 ^ 



In order to obtain estimate (5), let p be the radius of the ball B. Then 2p <G 
[M" 1 , M . ' : i for some ra. Now subdivide [M~ ,M^_ 1 ) in m n intervals each of 

length MJn. 1 - Lct 1 > bc sucn that 2p G [ZM,7\ Q + ^Kn 1 )- This means that 

TTi 2p ~ M " 1 < T 2p - 
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Using this estimate we know that each interval of J' n contains a ball of radius 
pj (I + 1), therefore a ball of radius p can intersect at most 1 + 1 + 2 such intervals. 
Using (6) we get: 



|l-Er, 



m < E \Jk-iS jl - in < E \ J J 

il)"' )3n 31)"' Jn 

< E M.^ 1 -^ <(l + 3)M-^ 

jl,— Jn 

If n is big enough, the fraction equals l e " + (1), which a priori may be unbounded; 
but we have: 



and: 

m: = cxp (logM„ 



log^ + n(logf -logC) 
logM„-logtf 



/logA/„-logtf + logtf (. q /, . „ 

= eXP ( logM^log/ ^ + n(log--lo S C 

^ CXP (( 1+ log^ g -log. ) + 

= exp f logM„-! f 1 l0g " V^^f-^)' 



logM„-logtf/ logM^.i 



and ejj — > as ri — -> oo. Then 

<m!i < (2p)- £ -. 

So we finally obtain p(B) < C|£?| 1_£ ™ _£ ™ . This estimate still depends on n therefore 
on but notice that e„ + e' n is monotone decreasing to 0, therefore if we fix n 
the inequality will hold for all B such that \B\ < M~ x . At this point it is easy to 
see that Ve > 3 6 such that any 5-ball B with S < 5 will satisfy inequality (5): 

p{B) < Cdiam(B) 1-£ . 

□ 

For all e > proposition III. 4 gives dim# J' > 1 — e, therefore 

(limn J > dim# J' = 1, 
that proves proposition III. 1 . □ 
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Now notice the very important fact that all estimates required in the proof are 
evaluated at points of the set J'; the assumptions of the statement can therefore 
be relaxed and we can require them to be true just on points belonging to such 
set. Moreover a corollary to the proof of proposition III. 1 is the somewhat weaker 
statement: 

Corollary III. 6. Let J n be as before. If mid > m^-d > 2 and rh~ n < Cm n , for all 
n then: 

J = O Jk dimnJ = 1. 

k£N 



In the following section we shall use only the weaker version we just stated. 



C. Reduction to the model system 

In this subsection we show that the results found in the previous subsection can 
be applied to the approximated bouncing ball system. We first build a set such that 
the dynamics of orbits that never leave this set is hyperbolic; in doing so we closely 
follow [Do07]. For small a we define the set 

B a = {(t,v) s.t. \4>(t)\ >a}, 

and compute the differential dF: 

dF=( 1 T » 

\2<f>{t + T{v)) l + 27»0(i + 7») 

Lemma III. 7. If v is big enough and F(t,v) € B a then dF is hyperbolic 

Proof. It suffices to check that Tr(dF) = 2(1 + T'(v)4>(t + T(v))) > 2, but for big 
enough v, T'(v) ^> 1, therefore, since by hypothesis \<p{t + T(v))\ > a we have 
hyperbolicity. □ 

Now we want to find an invariant cone field. In order to do so we find the direction 
corresponding to the expanding eigenvector in the limit T' — > oo, corresponding to 
high energies. We claim that a small cone around this direction is invariant for v 
large enough. In fact the eigenvectors of dF in the above limit are 

V+ = (St,6v) = (l,2$(t + T(v))) 
V- = {5t 7 Sv) = (1,0) 

Therefore, having fixed a small < c < a, the cone field defined on B a by the 
following expression: 



e t ,v = {(st,sv) s.t. 



< c} (7) 



is invariant for large enough v as V+ and V_ are well separated on B a . This means 
that if we take a curve whose tangent vectors lie in such cone field and we apply 
F we arc going to obtain (on the hyperbolic set) a curve whose tangent vectors 
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again lie in the cone field; moreover since the tangents to the curve are close to the 
expanding direction of the map, the dynamics will be expanding along the curve. 
For such reasons we now define 

r c = {(t,v) = (t,c + 24>(t))}, 

and curves from this family with big enough C will be good candidates for our 
purposes. 

The condition that an orbit never leaves the set B a controls the expansion rates of 
the maps F n . Now wc need to define a set to ensure the escaping condition. As we 
mentioned in the outline wc now define the following set: 

A e = {(t„, v n ) s.t. v n - u„_i > e} . 

As v n — tVi-i = 2</>(t„), we have that |20(t) > e X = A e . Wc can therefore select 
values of a, e and v such that there exists an interval J C S 1 such that 

J x {v >v} C A e D B a , 

and v is big enough for lemma III. 7 to hold true and for the cone field in (7) to be 
invariant for any v > v. Recall now the definition of F n : 

7(i) : S 1 3 t h-> (t, v = C + 2<p(t)) ir : (t, v) ^ t. 
F n = ir o F n o 7 . 

Lemma III. 8. For large enough C there exist positive constants Ci,C2,C_,C such 
that: 



ViG f] F^(Jo), ^HiC + ek)^ 1 < \F^(t)\ <C 2 \{ (C + 3Ak) 



7-1 



fe=0 



k=0 



k=0 



Proof. Let F k (t,j(t)) 4= (tfc,i>fc). First we get an estimate for d ^ 1 ; we know that 
tk+i =t k + T(v k ), then: 

— — = 1 + T K) — . 

dtk dtk 



If t S P| F k (J) and C is big enough, we know that the cone field C is invariant i.e. 
(1, ^) G 6, therefore 



that implies 



dtfc+i 
dt k 



dt k 



mtk+i) 



< c, 



(1) + Const wjp^tfc+l) + (c)). 



and since (t k , v k ) € B a and c < a wc can find positive C\ and Ci such that 



dtk+i 
dt k 
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Now since each Vk) € A e we have bounds on v^, i.e. 

(7 + ek < v k <C + 3Ak, 

where C_ and C are respectively the minimum and the maximum v of the curve T 
and by (1), A is the maximum of </>(t). Therefore 

< c 2 (c + a^fc) 7 ^ 1 . 

Then by the chain rule we get the required inequality. □ 

In order to apply the results found in the previous section we have to ensure that 
CiCp- 1 • | J 1 > 2; this can always be obtained by choosing C large enough. Now 
we can apply corollary III. 6 and obtain dimfl-(<? fl J) = 1. □ 

As a final remark we would like to point out that the same idea works for the non- 
approximated map with some slight modifications; the proof has not been included 
in this paper for sake of simplicity. 



Ci (C + ek) 



7-1 



< 



dt 



fc+i 



df* 



IV. EXISTENCE OF ELLIPTIC ISLANDS FOR ARBITRARILY HIGH 

ENERGY 

In this section we are going to prove Theorem B: 

Theorem B. Let a < 1. Then in the static wall approximation, for almost all 
sinusoidal motions there are elliptic islands of fixed period 2 for arbitrarily high 
energies. Moreover if2/3<a<l we can prove the same result for all sinusoidal 
motions. 

First we recall the definition of elliptic island: if an elliptic fixed point p for a 
two-dimensional symplectic map F is surrounded by a invariant set of closed curves 
on which the dynamics is conjugated to an irrational rotation on the circle we say 
that p is surrounded by an elliptic island. Such islands are obviously Lyapunov 
stable. 

The outline of the proof of Theorem B is as follows. In section IV A and IV B 
we build a reversor map by exploiting a symmetry of the system; we recall that 
a reversor is an idempotent map that conjugates the dynamics with its inverse. 
Following a standard technique in the theory of reversible maps (see e.g. [LR97]), 
we use the locus of fixed points of the reversor map to find a number of periodic 
orbits; most of them will be hyperbolic but by fine-tuning the amplitude A of the 
oscillation of the plate we can turn some of them into elliptic periodic points. In 
the main approximation (defined in section II) it is easy to state the cllipticity 
condition (section IV C) in terms of A. We can actually find conditions to ensure 
that the multiplier of such periodic points belong to some given sub- interval of S 1 . 
This fact will turn out to be useful when dealing with resonances. Such conditions, 
along with a non-degeneracy condition on the Birkhoff normal form that we prove 
in section IV F, is sufficient to establish the presence of an elliptic island around the 
periodic point (see for instance [La93] or [dlL99]). The cllipticity condition (section 
IV D) turns out to be an arithmetic condition on the parameter A that is satisfied 
(section IV E) by infinitely many periodic points for a set of full measure of A for 
all weak potentials. The same proof gives the stronger result that for 2/3 < a < 1 
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the statement is true for all parameters A. 

The techniques developed in the following sections allow us obtain some interesting 
quantitative results. As an example, in proposition IV. 15, we state conditions on 
7 guaranteeing that the Lebesgue measure of the elliptic islands obtained with our 
construction is cither infinite or finite. 



A. Construction of periodic orbits 



Recall first the definition of F: 

zr f 1 t + T(v) 

F: \v " v + 2Ht + T{v)) T (v) = C-v^. 

We are going to find periodic orbits by building a set with special dynamical 
properties and then considering intersections with its forward and backward images. 



Remark 1. If <p is odd with respect to some point to, then F has a reversor map 
R such that 

i? 2 = Id RFR=F~ 1 . 

We can explicitly write R as follows: 

R: (t,v) i * (2t Q -t-T(v),v). 

If we square the map we get the identity and it is an easy check that R conjugates 
F with its inverse. Notice also that being defined on a cylinder, if <j> is odd with 
respect to to it has to be odd with respect to to + h as well. 

We are going to define the set £ (for locus) of fixed points of R: 

e = {(t,v)s.t. R(t,v) = (t,v)}. 

The set I is the disjoint union of the following two curves: 




Such curves wind around the cylinder as v increases. It is more convenient to 
partition such curves in pieces that wind just once around the cylinder in order to 
get graphs of (single valued) functions of t. This can be easily done by inverting 
the 1-1 map v <—> T; let this inverse be v(T). Define now: 

VneN £ n (t) = (t,v(2(t -t)+n)). 

The curve t+ corresponds to even values of n while to odd values. Subscripts 
will always refer to branches and superscripts will always refer to iterates of the set, 
i.e. for keZ,£l = F k l n . 

The important dynamical property of I is that 

Vxe£ F k x = RF~ k Rx = RF~ k x, 

therefore, if F k x belongs to I as well, we get F k x = F~ k x that implies that the 
orbit of x is periodic of (possibly not least) period 2k. Therefore points belonging 
to l k n I for k ^ are periodic points. The issue is now to figure out when the 
corresponding periodic orbits are elliptic or hyperbolic. Taking inspiration from 
[GLOO] we work out from scratch the period 2 case. 
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B. Period 2 orbits 

First we classify period 2 orbits. This turns out to be quite simple, as the following 
proposition shows. To fix notations, let {pi,p 2 } be a 2-periodic orbit, p\ = 
and p 2 — (£2,^2); since v = v 2 , we have 0(ii) = — 0(t 2 ). 

Definition IV. 1. Being is a sinusoid, 2-periodic orbits can only be of one of the 
following types: 

• 4>{ti) = 0(^2), such an orbit will be called a (+)- orbit; 

• 4>{t\) = — </>(i 2 ), such an orbit will be called a (-)-orbit; 
Proposition IV. 2. Let {pi,p 2 } be a 2-periodic orbit; there can be two cases: 

• _Pi>P2 G the orbit is a (+)-orbit; 

• T(vi) = T(v-2) = h mod 1; the orbit is a (-)-orbit. 

Proof. Let us write the condition for p\ = (t\, vi) and P2 = (£2, ^2) to be periodic: 

h 1 1 ti+T(yi) = t 2 

t 2 1 ^ ii +T(ui) +T(u 2 ) = *i 

«i h-> ui + 2^(t + T(vi)) = v 2 

v 2 ^ «i + 2 (<£(i + T(«i)) + </>(t + T(«i) + T(«a))) = «i. 
So that we have the two conditions: 

T(«i) + T(«a) = modi 0(t + + = 0. 

As 0(f) is a co-sinusoid, the second condition can be true only under one of the two 
following assumptions: 

• ti + T(vi) = 2to — t\ therefore p\ £ £; the same is true for P2 

• h + 7>i) = h + i therefore T( Vl ) = T(v 2 ) = \- 

Notice moreover that as <j> is odd with respect to to, <j> is even with respect to the 
same point to, therefore orbits of the first kind are (+)-orbits. On the other hand, 
orbits of the second type satisfy the opposite condition <f>{t\) = —(f>{t\ + 1/2 = t 2 ) 
and so they are (— )-orbits. □ 

C. Elliptic locus for period 2 orbits 

In this section we find a subset of the phase space such that all 2-orbits that lie 
in the set are elliptic. 

Proposition IV. 3. Let the two points of the orbit be (ti,v\) and (t 2 , v 2 ); let T[ = 
T'{vi), but 4>\ = 4>{t 2 ), 4>2 =f <P(ti) and similarly for (f>. Notice that by proposition 
IV. 2 we have 4>i = ±</>2- Let i>i = (\>iT[ and —1 < c\ < c 2 < 0; define the following 
sets: 

E Xc 2 {(V1,V2) S.t. V1+V2+ V\V2 G (Cl, C 2 )} 

E c lC2 {(^1,^2) S.t. V\V 2 G (c 2 ,ci)} 

Then (±) -orbits belonging to E^ lC are elliptic with multiplier A such that 3?A G 
(l + 2ci,l + 2c 2 ). 
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Proof. We start by computing the differentials dF on each point on the orbit: 



dfi = 



1 



24>x l + 2<t> x T[ 



dF 2 = 



1 



J 2 



202 l + 2(/) 2 T^ 



The condition for any matrix M to be elliptic is that |TrM| < 2; moreover, for any 
elliptic matrix M . ^TrM is the real part of its multiplier. By direct computation 

of the trace of the product of the two differentials and using the (±) relations on <f> } 
we obtain: 



-Tr (diW 2 ) = 1 + 2 (faTl + foT'i 



2T[T 2 <f>\(t>2- 



(+) 



-Tr (diWa) = 1 + 2T 1 / r 2 '0 1 2 . 



(-) 



By direct computation we obtain the following conditions in i/j that ensure cllipticity 
and the supplementary condition on the multiplier: 



V1+V2 + V\V2 G (ci, ox) 



(+) 



V\V2 G (ci,c 2 ), 

that are the defining conditions for the sets E^ c . 
Notice that since 4>i = ±</> 2 we have that 



Since V = T 1 - 1 /7 an d 



1/2 



jX/t _ j>Vt 



rpl ■ 



A = max 



we get 



(-) 
□ 



v 2 



± 



±1" as Ti 



Even if the shape of the sets E~^ c is not very complicated, it is convenient to state 
a sufficient condition in terms of just one parameter v. Fix e small, then if we let 
vi < V2 big enough, we have (1 — e)|^2| < \vi\ < \v2\] a direct calculation yields 
the following sufficient conditions for 2-orbits to be elliptic satisfying the required 
condition on the multiplier: 



(+)-orbits: v 2 G (ci,c 2 ) (-)-orbits: \v 2 \ G (|c"|, \c 2 1) 



(9) 



where c'^c^ and c'{,c 2 are e-close to ci and C2- Note the we can get all elliptic 
orbits just by taking ci = — 1 and c 2 = 0. 



D. Description of I and ellipticity condition 



In this and the subsequent sections all pictures and geometric constructions are 
made keeping in mind the coordinates (t, T(v)). In such coordinates I is represented 
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by a straight line, and it is much easier to have geometric intuition about the 
dynamics. Recall that the simple choice for <fi given by equation (1) implies the 
following expression for <fi: 

4>{t) = Acos(2irt). 

Now we claim the following 

Proposition IV. 4. Let v(T) be the inverse function ofT(v). Then there exist real 
positive numbers C\ < G% such that, for any n <m€N the following condition 

( 1 C 2 1 \ / 1\ „ „ / 1 Ci 1 \ 
v m- r - rr <v[n+-)+2A<v[m-- _ 

implies the existence of a 2-periodic elliptic point close to T = n +1/2 and T = 
m + 1/2 such that its multiplier satisfies SRA € (1 + 2c\, 1 + 2c-i). 

Proof. Consider l\: let Ufc(r) = v (2 (r — to) + k) so that: 

4(r) = (r,i;„(r)+20(r)). 
So for A = this is just a line in the (t, T(u))-planc. As A increases, such line 




Figure 2: On the left the reference picture for A = 0, on the right the situation for A > 

deforms and resembles the shape of <f>, as we can observe in figure 2. As we noticed 
in proposition IV. 2, 2-orbits obtained by intersecting £ and I 1 are (+)-orbits; we 
now claim that the highest energy point ^2,^2) of the orbit lies where <fi > 0. In 
fact we know that 4>{t% ) is the velocity gained on the collision at time t\; since 
we want that V2 > v± we need such quantity to be positive. Having that fixed, 
vi = 4>(ti)T' (V2) = <f>{t2)T' {V2) and v\ < V2, therefore we have that condition (9) is 
satisfied if (t2, 1*2) £ ^ H I 1 belongs to this set: 

First notice that this set is an (l/(A ■ T'(w)))-thin strip that lies {I /{A ■ T'(v))) 
on the right of the vertical line r = (that corresponds to c' 2 = —1). By direct 
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inspection we obtain that in (t, T) coordinates, each branch of £° is a straight line 
with fixed angular coefficient -2 and each branch of £ near r = is approximated by 
a parabola that intersects r = with positive derivative (close to 2) ; the maximum 
of such parabola is given by the equation: 

0(Tj ~ T'{v n (r)) < T'{v) < T'{vY 
Figure 3 illustrates the properties we just described. The key fact to notice is that 



T = 




Figure 3: Explicit construction, in (t, T) coordinates, of values of the parameter A for 
which we have an elliptic periodic point of period 2 with given bounds on the multiplier. 

the values of A we arc seeking arc close to values A of the parameter for which the 
intersection lies on the vertical line r = 0. Let us compute the intersection of £ m 
and t\ with the vertical line r = 0: 

£ m (0) = (0,«(-l/2 + m)) (10) 



£(0) = (0,«(l/2 + n)+2A) (11) 

Therefore if we want £ and I 1 to intersect on the line r = we need to find A such 
that points in (10) and (11) are equal for some n,m, i.e.: 

v(n+±)+2A = v(m-^) (12) 

Now it is clear that we can find A\ and A% as in the picture such that the intersection 
lies on the boundary of c , . Using the properties we described above it is also 
clear that, in (t, T) coordinates, the distances between intercepts of £ Y corresponding 
to each A\ with the vertical r = are linear functions of the t coordinates of 
the intersections themselves, therefore of order Q (I / (A ■ T' (v))) . More precisely, 
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mimicking equation (11) and recalling that T" ~ T 1 1 / 7 , we obtain that there exist 
C\ and C2 such that if 

f 1 C a 1 \ ( 1\ „„ / 1 C x 1 \ 

» m-- ; 1 : ; <u n+- +2 J 4<w m — - . . 

\ 2 A m x ~ x h J \ 2 J \ 2 A m 1 - 1 h J 

then the intersection (£ m nfi)nM,^8. □ 

The condition we obtained is basically an arithmetic condition on A and 7. In 
the next subsection we prove that this condition is satisfied for parameters A as in 
the statement of Theorem B. At that point we will be only left with checking the 
non-degeneracy condition. 



E. Arithmetic condition 



In this section we are going to prove an arithmetic result that is of independent 
interest; for simplicity we state the arithmetic condition in a slightly simplified form 
with respect to the case in consideration. Namely we drop the 1/2 that appears in 
the statement of proposition IV. 4 and we reverse the signs of Cx and C2. One can 
easily verify that this does not affect the proof in any sense. 

The condition is reminiscent of the Khinchin's theorem on Diophantine approxi- 
mation [Kh64]. In fact we want to investigate parameters 7 and a such that the 
following inclusion is true for infinitely many n and m G N: 



2 -f 
—m 5 



for an appropriate (and fixed) choice of C2 > C\ > and £ > 0. In our case a = 2A 
and £ = 1 — I/7. 

Let us first introduce some useful definitions: 

Definition IV.5. Let us fix £ > 0, 7 > 1, C 2 > Cy > 0. Then 



C 



' , m 



2 -t 



'a — II ^a,m & a ~ <j a 3 



met 



X 



a = [n lh + a, n G n} . 



Using this notation a parameter a satisfies the arithmetic condition if the cardi- 



nality 



x a n sr, 



is infinite. 



Definition IV.6. Let n, k G N: 



a G R s.t. 



} 



# = {s£l+ s.t. (n lh + a y G %\ 



= {«el + s.t. (n 1 /"' + a y g %, (n+k) } ; 
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Clearly si n = \Jk^k> moreover if si n « = \J n>na #i n , then si = f) si n ° = 
limsup n ^ 00 ^™ 

Lemma IV. 7. For all real 7 > 1 si is a residual set in R + . 

Proof. Each si n ° is open since it is a union of open sets. Moreover it is dense since 
the distance between endpoints of consecutive intervals belonging to Sf a goes to 
as m — * 00, and so do the distances between endpoints of the intervals belonging 
to si n as n — * 00. As the point is the limit point of the left endpoints of the first 
interval in si n we conclude that s/ n ° is dense in M + . □ 

Let us define the following conditions involving £ and 7: 

f < 1 (diverging); 



£ < — (overlapping). 

7 

Notice that since 7 > 1 the overlapping condition implies the diverging condition. 
Now we can state the result as follows: 

Theorem IV. 8. If the diverging condition does not hold, then si has measure 0. 
If the diverging condition holds then si has full measure in R + ; moreover if the 
overlapping condition holds as well then si is the whole R + . 

Notice that, as in Khinchin's theorem, we obtain that the required property is 
satisfied either by a null set or by a full measure set. This dichotomy seems to be 
quite common in approximation problems similar to the one we arc studying. 

Proof. The proof will be presented in four steps. 

a. Step one By restricting to a compact interval of possible parameters we get 
rid of the parameter a in the definition of the sets Sf 

We restrict ourselves to a bounded interval of parameters a £ [a, [3] with the con- 
dition 0/a < C2/C1. This condition implies the middle inequality in the following 
expression: 

l Oi Oi L/2 O2 O2 

(3 a a [3 a a 

We can therefore build a superset and a subset ^ a p* such that £f *g D 5f a Z> ^ a /3* 

for any a £ [a, [3], therefore proving the result using either Sf*^ (to obtain estimates 

from above) or Sf a/ 3* (to obtain estimates from below) provides it for the original 
problem as well. Therefore we now fix once and for all < a < [3 < (C2/C\)ol and 
we prove the result for parameters a in such interval; to avoid the introduction of 
new notation we redefine si, si n , si^ 1 as their intersection with the interval [a,/?]. 
As we can write R + as a countable union of such intervals we obtain the desired 
result. 

b. Step two We investigate the structure of the sets si n 

Define <5£, A£, AJJ, I£ as in figure 4. The following lemma provides some useful 
estimates: 
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a 



< +1 k 



of" a/" 



JL 



Figure 4: Definition of tfjf, A£, A£, 



Lemma IV. 9. Define the following positive quantities: 



C*2 Ci 



. C*2 Ci 



a P (3 a 1 

and the following set K n = {ieN s.t 7^ 0}. Then: 

81 = -(n + h)-*--*- 1 - 1 ^ + h.o.t. forle U + ,t 
7 



(El) 



A£ = -(n + k)-< 1 - 1 W + h.o.t 



(E2) 



= -(l--\kn 1 ^- 2 + h.o.t. 



(E3) 



K n ~ [7a • n 1 " 1 / 7 + h.o.t.,j(3 ■ n Y - lh + h.o.t] n N 
Proof. We first bound the length of the intervals s4£: 

81 < (n+k + ^-in + k)-^ H - (n + k+^-(n + k)-t 



1/7 



— (n + k) 

7 



-«-(l-l/7) 



h.o.t. 



(E4) 



The bound from below is similar and yields the expected result. Then we estimate 
the length of the intervals 1%: 



A£ = (n + fc + 1) 



l h _ „i/7 _ 



(n + k) 



l h _ „l/7 



(n + /c + l) 1/7 - {n + k) lh 

i(n + fe) 1/7_1 + h.o.t. 
7 
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Next the offset of two subsequent srf£ : 



a: 



(n + k + 1) 1/7 - (n + 1) V7 - ( (n + fc) V7 - n 
(71 + fc + 1) 1/7 - (n + k) lh - ((n + I ) 1/7 - n 1 / 7 ) 



1/7 



\ x h _„ih 



(n + k) 



• W7-1 _ n V7-l 



h.o.t. 



.1 (l--) k-n 1 ^- 2 + h.o.t. 
7 V 7/ 



Finally we estimate K, 
K n = 

notice that 



n lh + a 



7 ( , , \ 7 

1/7 



n, n i/7 + /3 



nN 



1 ° 



■ n 1 - 1 /" 1 + h.o.t., 7/3 • n 1 " 1 h + h.o.t. 



771 



1-1/7 . 



(ft -a). 



□ 



c. 5iep £/iree Overlapping regime 
From the previous estimates we can already obtain the result in the overlapping 
regime. In fact (E4) implies that k is 6'{n 1 ^ 1 ^'), therefore, by (E3), A£ is 
G{n- X ). This implies that if -£ - (1 - 1/7) > -1 (i.e. £ < 1/7 that is the 
overlapping condition) sd£ and ^/ fe ™ +1 will eventually overlap as, by (El), the 
length of the intervals si£ goes to zero slower than their offset. Since they overlap 
and they are moving like 1/n, they will eventually get out to the left so for each 
fixed k they are going to cover the whole interval, therefore sf n is going to cover 
[a,/3] infinitely many times, and s4 will contain [a, /?]. 



d. Step four Non- overlapping regime 
We will now focus on the strictly non-overlapping regime i.e. £ > I/7; the critical 
case £ = I/7 will be considered later as it is just a combination of other cases. 
Define, for any Borel set E C [a, 0], F(E) = J£(E)/((3—a) (the normalized Lebcsguc 
measure) as a probability measure on [cc, /?]. We are going to prove that the set 
si has either full measure or measure zero using the following strong form of the 
Borel-Cantelli lemma: 

Lemma IV. 10 (Borel-Cantclli-Erdos-Renyi [ER59]). Let {^4} be a sequence of 
events on a probability space (f2, <^*,P). If 

00 

V^P(^,) < 00 (convergence), 
n=l 

then P(limsup.e/„) = 0. If instead 

00 

P(«g^) = 00 (divergence) 

71=1 
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liminf /T _^ n ^3- = (weak independence), 



En 
k,i=i r \^km 

(ELiPK)) 

then P(limsup s e/„) = 1. 



We will get the result by verifying the above conditions for the sets si n . The fol- 
lowing two lemmas deal with the estimates of the convergence/divergence condition 
and the weak independence condition respectively. 

Lemma IV.ll. Fix n and let N tend to infinity. Then: 

N 1 

Y^P(^") = - — -IN 1 -* + h.o.t. (13) 

n—n 

Proof. Let us first compute: 

JSf(^) = J2 E = ~ E (n + fc)- £ - 1+1/7 

k&K n k&K n ' kGK„ 



n~t- 1+lh + h.o.t. = ([3 - a)£n- 6 + h.o.t. 



7 

' k£K 



Then by normalizing and summing on n we get (13): 

JV N 



□ 



Observe that the diverging condition implies that (13) diverges as N — > 00, while 
the series converges if the condition is not satisfied. 

Lemma IV. 12. Fix n, let N go to infinity, then: 

N N 

V V ¥(^ n £/ m ) = e 2 -, -n 2 - 2 ^ + h.o.t. 

^ ^ v ' (1 - £) 2 



n—n m—n 



Proof. First of all, by symmetry we can consider just m > n by paying a factor 
of 2 and some diagonal higher order terms; then we separate again srf n in their 
respective components srf£ 



N N N N 



^2 E p K n ^ m ) = 2 E E p K"^ m ) + k°.t. 



n—n m—n n—n m—n 



N N 

= 2 E E E E p K<) 

n=n k£K n m=n lGK m 
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Now let us focus on the last two sums, fixing momentarily n and k. We can write 
the following equality: 



N 



Y p (^k^r) = J2 P (■ s/ k B p k ) + err ° r term 



(14) 



m—n[^K m p— 1 

where we define the sets B^ k as the p-th wave (see figure 5): 



B': 



P£ is an appropriate number that is estimated by the next lemma and the error 
term is due to the fact that the last wave could be incomplete. 



Lemma IV. 13. 



The following estimate holds 



r>n 



i •fc + om 



4+p and is A£ +p . 



is 



Proof. As it follows from (E3), the offset between 

the union over all r such that s^^+p intersects 1%; for each fixed p there is a wave of 
.e4 r +p that will be inside 1% for some time and then leave the set. PJ} is the number 
of waves that will pass through Tfc in the time N. This means that: 



(N + k + P 7 h !) lh - N lh ~ (n + fc) 1/7 ~ « 1/7 , 



that is: 



iV 1 ^-! (fc + P fc n ) - n 1 ^ 1 /,- => P fe " = 



0(1) 



□ 



The error term in (14) can be easily bounded by the order of as we 

miscount of at most 1 wave. We are left with the computation of P {srf£B™ k ) . Each 
B£ k is the union of intervals that are A£ + apart and S^+p l° n g- The ratio of such 
data gives the portion of the interval covered by each wave. 

Lemma IV. 14. Let us introduce the parameter r\ = (k + p)/k. Then we have: 



k+p 



A 



k+p 



1 k 



rj -i- 1 



+ h.o.t. = AjJ(p). 



a 



Bf 



"k+1 



k+p 



Figure 5: Relation between waves Bp and sets srf£ 
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Proof. By definitions of 5^ and A£ we have: 

_ g 7 (r + k + p)-Z-^ 
Al +p 7-1 (fc + p)(rV7-2) • 

We need an estimate on r: acting as before for the computation of P wc get: 

•y 

r 1 /i r - 1 (fc+p)~n 1 /T- 1 * r = n('^) 7 " 1 . 

We rewrite the previous expression as: 

7 (nrf'^-^ +r 1 k)-Z- 1+1 /-> 



J-lTjk- nV7-2 . r/ (7/(7-l))(l/7-2) ' 

that is: 

(( n + ?r l/(7-l) fc ) r? 7/(7-l))-«- 1 + 1 ^ „l-£ 

£ — — — — - — LI — 7r ^ + h.o.t. = l—L- -—77— + h.o.t. 

7-1 k ■ n x /7-2 . ?r 7/(7-i) 7 - 1 fc 

□ 

Therefore P (&/£B£ k ) = AJ?0) • P(^) + ^(A£ +p /(/3 - a)), where the error term 
comes from the granularity of the set B. Therefore: 

pn p n 

E F W) = p«)E^(p)+ii.o.t. 

p—1 p—1 

pn 1 — f-y 



1— f fc 

= PQQl 7 1 (1 — E( fc +P)~+ W - 

~ A; 7-1 p=1 

1 * / / (Mh\ (l-£k\ 

= P«")£ l (i _ -(((fc + P») — )-fc— ) + h.o.t. 
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H^) h , ( (-13177) 
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fc ^„W.lW-lj + h.o.t. 
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Now we sum over n and k: 



N 1 //v 1 -? \ 

E E Wrn» H — - 1 ) + h -°-t- 



AT 



n—n 



k£K n 



1 w 



£ 2 N 2 ~ 2 ^ + hot 

2(1 -0 2 + n -°- 1 - 

Recalling the factor 2 we had in the beginning of the estimate we get the desired 
result. □ 

The last two lemmas prove that the weak independence condition is always sat- 
isfied independently of £. Therefore we have only to check the diverging condition. 
In the diverging regime we can conclude that the set srf has full measure, as we 
found a subset of full measure; in the non-diverging regime we can as well conclude 
that stf has zero measure, as it is contained in a zero-measure set. 

e. Step five The critical case £ = 1/7 
For £ = I/7 we have that the overlapping condition is satisfied for small enough 
a, as 5% grows bigger as a decreases. As we notice from lemma IV. 14 we can find 
a critical a such that for a < a we have overlapping and for a > a we have no 
overlapping. □ 

In the case in our consideration £ = 1 — 1/7; this implies that the diverging 
condition is always satisfied and so we proved Theorem B up to the non-degeneracy 
condition. 

Remark 2. For this particular value of £ the overlapping condition reads 7 < 2; 
the critical case is therefore 7 = 2. In such case it is actually possible to explicitly 
compute a value A £ whereas in general for 7 > 2 we still ignore if M. + \ si ^ 0. 

Remark 3. The technique we developed can be applied to (— )-orbits as well. The 
arithmetic condition relative to such orbits turns out to be more restrictive than 
the one for (+)-orbits; more precisely we get £ = 2 — 2/7. This implies that the 
diverging condition is not anymore guaranteed. In fact it fails for 7 > 2, which 
means that such orbits appear for arbitrarily high energies for almost all sinusoidal 
motions only for 7 < 2. 

Having studied all possible 2-periodic orbits, we notice how the conditions we stated 
are actually also necessary conditions for the presence of elliptic 2-periodic orbits. 
This implies the following interesting results: 

• by the previous remark, for 7 = 2 there is at least a value of A such that the 
system has only finitely many elliptic islands of period 2. 

• if 7 > 2 we have infinitely many (— )-elliptic islands only for a null-measure 
set. Notice however that lemma IV. 7 does not depend on £, therefore such 
set is non-empty. 
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Remark 3 allows also to prove the following proposition regarding the total mea- 
sure of elliptic islands, given the following estimate of the measure of a single elliptic 
island: the measure of each (-f)-elliptic island is of order 1/T' 3 . We are going to 
prove such estimate in proposition IV. 18. 

Proposition IV. 15. If 7 > 4/3 the total Lebesgue measure of (+) -elliptic islands 
of period 2 is finite. If 7 < 4/3 the total Lebesgue measure of elliptic islands is 
infinite. 

Proof. We can obtain a rough upper bound to the total measure of elliptic island 
by summing the measure of a single island over all possible islands that can appear. 
As there could be one for each T ~ n + 1/2 we have the following estimate: 

00 00 
Leb(islands) < Const • ^ ^ = Const ' Yl ™~ 3(1 ~ 1/7) - 

n—l n n—1 

The series converges for 7 > 3/2. To sharpen the estimate we need to take into 
account that for some of the n we do not have an elliptic island. For 7 < 2 this 
may be estimated quite easily, as we have the following expansion: 

(ny + A)i = n + ^An 1 - 1 ^' + o(l). 
From the previous expansion it is clear that : 

{(r^ +^1) 7 } = {jAn 1 -^} +o(l). 

such function has an infinite number of branches, let us denote such branches by k. 
Each branch will start at ~ fc 7 ^ 7 ^ 1 ). The arithmetic condition can be expressed 
in terms of k in the following way: 

{n 1 ^} <0(fc- 1 ). 

Figure 6 pictures this condition; for each branch k we have elliptic islands until 
the fractional part grows too large and the arithmetic condition no longer holds 
true. Given this fact it is easy to estimate the number of islands belonging to the 

1 — 1. 

/c-th branch. We compute the derivative of n t for rik+i, obtaining a linear lower 

l— — 

bound on the growth of n 1 in the /c-th branch: 

> • (n - n k ) ~ (k + . ( n - n k ). 

The smallest n for which the arithmetic condition fails can therefore be bound from 
above by requiring: 

(k + l) - ^ ■ (n - rife) < (k + i.e. n - n k < (k + . 

therefore for the fc-th branch we get at most ^/c^ 1 ~\ elliptic islands. Wc now 

multiply this number by the measure of such islands and sum over all branches k 
to get the total measure: 

00 

Lcb(islands) < C ■ ^ k^k~ 3 , 
fe=i 

that converges for 7 > 4/3. Notice that along the same lines we can get a lower 
bound of the same order, that means that the total measure of islands diverge for 
7 < 4/3. □ 
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F. Non-degeneracy condition 

According to general KAM theory, there exists a stability island around each 
point of a periodic orbit provided that generic non-resonance and non-degeneracy 
conditions are satisfied. Following [La93]: 

Definition IV. 16. An elliptic fixed point p of a two-dimensional symplectic dif- 
feomorphism / is said to be general elliptic if: 

• the multiplier X p is such that A£ ^ 1 for k = 1,2,3,4 (non-resonance up to 
order 4); 

• the Birkhoff normal form is non-degenerate, i.e. a quantity that can be written 
in terms of derivatives up to fourth order is different from zero (see below). 

Theorem IV. 17 (KAM). If p is general elliptic, then it is stable, i.e. for each 
neighbourhood U of p there exist another neighbourhood V such that Vfc, F k (V) C 
U. 

Stability around the point implies the presence of an elliptic island. As noted 
before, the construction we described yields elliptic points with multiplier which 
can be chosen to belong to some prescribed interval; this implies that we can a 
priori avoid resonances. The non-degeneracy condition can be explicitly computed 
by the following procedure: 

• We perform a linear change of coordinates such that the differential of the 
map dF at the fixed point is a rotation in the new coordinate (£, rf) where 
(£ = 0,77 = 0)^ p. 

• We compute the Taylor expansion coefficients up to order four (excluded) in 
the coordinates u = £ + in and u obtaining the following expression: 

u i ► \ p u + A 3 u 2 + A±uu + A 5 u 2 + A 6 u 3 + A 7 u 2 u + A 8 uu 2 + A 9 u 3 + (4) . 
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Figure 6: The arithmetic condition for 7 < 2 
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We compute the following expression: 



3 



i^X p A 7 )+3\A 3 \ 2 ^ + \A 5 \^ 



The non-degeneracy condition requires that to ^ 0. The coefficients Ai contain 
derivatives of T up to order 3 and derivatives of <f> up to order 4. As for high 
energies we have T[ 3> 1, instead of computing all Ai exactly, we perform an 
expansion in terms of powers of T[ and compute the highest order non-zero term, 
taking into account that ellipticity implies </>T/ = Vi ~ 1 (i.e. condition (9)). We 
find by direct computation[l] that the highest nonzero term in uj is of order T[ . As 
a further simplification we note that we have T 2 = T[ + (T"); this implies that if 
we compute uj by setting T[ = T' 2 = T' and find a quantity bounded away from zero 
in this limit, it will be bounded away from zero also for all T[ sufficiently large. The 
coefficient of order 3 turns out to be the following polynomial in v = <f)T' e (—1,0): 

1{v 2 + Au + 6) + \-iv{2 + vf + |^(3 + v) 2 ]. J . 

where 

D = ^2(l-Rc(X p -d z (F 2 (p)))). 

is of order \[T' and d z (F 2 (p)) is the holomorphic derivative of F 2 with respect to 
z = t + iv. Notice that the fractions involving the multiplier A p = cxp (i8p) give 
respectively i ■ cot(0 p /2) and i ■ cot(30 p /2). Of course 9p is not independent of v, 
but recall that since we can control the multiplier, we can assume both cotangents 
to be bounded away from zero and positive. It is easy to check that for v £ (—1, 0) 
such polynomial is bounded away from zero, as each term in the sum is positive. 
This is enough to establish the presence of an elliptic island around each periodic 
point found by the construction, proving Theorem B. Notice that the expression for 
a>3 does not involve derivatives of ip of order higher than 3 and derivatives of order 2 
and higher of T as such terms appear only in terms of lower order in the expansion 
in T' (see below). To conclude this section we now prove an estimate regarding the 
size of the elliptic islands we obtained. 

Proposition IV. 18. Consider a 2-periodic orbit of type (+), given by the points 
(ti,t>i) and (£2,^2) and such that the multiplier is bounded away from resonances of 
order up to four; we define T' = (T'(vi) + T'^)) /2. Elliptic islands of type (+) 
around such points have area of order (T') -3 for big enough v. 

Proof. We consider the map F 2 expressed in terms of the variables u, u defined 
above, close to a periodic point p; for simplicity we assume p = 0. Recall that the 
variables u and u are related to z = t + iv and z by a linear symplcctic transfor- 
mation, i.e. z = biu + 62U and bibi — 6262 = 1. In such variables one can write the 
map as follows: 

A 1 (u,u)u + A 2 (u,u)u where Al(0,0) = A p , A 2 {0,0) = 0. 

One can obtain all terms of the Taylor polynomial of F 2 in such variables by appro- 
priately differentiating the functions Ai and A 2 with respect to u and u. We claim 



{2 + v)4> 

0J3 - 



64D 6 /T 
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that the term of order n is of order at most (T') 3+ ( n *)/ 2 . By direct computation 
we find that 

A 3 , M, A 5 ~ (T'f 2 A 6 ,A 7 , A 8 ,A 9 ~ (T') 3 . 

This, along with the estimate we claim, is sufficient to prove that the area of the 
elliptic island is of the required order. In fact, one can perform a rescaling u i— > Au, 
obtaining the following (symbolic) expression: 

oo 

u i-> X p u + J2 4 (n) A n_1 {"> «}"• 

n=2 

Therefore by choosing A such that A^A n ~ 1 < 1 we obtain that the linearized 
part is dominant in a disk of radius of order A around the origin. The explicit 
computations and the claim allows us to take A ~ (T')~ 3 / 2 ; the result follows by 
recalling that the map z i— > u is symplcctic and therefore it preserves the area form. 
We are now left with the proof of the claim, i.e. to prove that A^ n > < (T") 3+ (™ -1 )/ 2 . 
First we obtain by direct computation a relation between the coefficients a% of the 
Taylor expansion in terms of z, z and the coefficients Ai of the expansion in terms 
of u, u 

z i— > di{z, z)z + a 2 (z, z)z z — biu + b- 2 u. 



A\ = bibidi + bib 2 a 2 - b\b 2 a 2 — b 2 b 2 a\ A 2 = bib 2 a x + b\b\a 2 — b 2 b 2 a 2 — b\b 2 a\. 

We can obtain all coefficients A^ n ' by applying the relative differential operator to 
the appropriate Ai\ the key fact to notice is that, bi being constant, the differential 
operator will operate only on the a%. One can express d u and <9 fi in terms of dt and 
d v in the following way: 

a h+b 2 b 2 ~bi bi + b 2 ~ b 2 -bi 

Ou = — - — o t + i — - — d v d u = — - — d t - i — - — d v . 

We are going to explicitly compute the coefficients of dt and d v to check that they 
are respectively of order T'' 1 / 2 ) and T'~( x / 2 \ Then we find a general expression 
for the order of arbitrary derivatives of ai and a 2 . Explicit calculations provide the 
following values for af 

01 = f 1 + 2<j>{T[ + T 2 ) + 2<f> 2 T[T 2 \ + i (2$ + 2<j) 2 T[ - TL^iTk _ <$T[Tl 



a 2 = ( -2ij)T 2 - 2'(j) 2 T[T 2 ) + i ( 20 + 20 2 T 1 / + ^4~^ + 4>T[T! 



2 

Where recall that we defined T[ — T'(vi); by defining as before 



D = ^2(l-»(Ap-oi(0)) ~ T'^' 2 \ 
and T = {l/2)(T{ + T 2 ) we compute 61 and b 2 : 

h = (2(f> + 24> 2 T[ — T' — ii>T[T 2 + 3A p ^) 

b 2 = I (-24, - 24> 2 T{ — T' — 4>T[T 2 - i (2^ + 24?T{Tty . 
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From the above expression, recalling that <f>T' ~ 1, we notice that each hi is of order 
T'^ 1 / 2 ). Now wc compute the coefficients that appear in the change of variables 
(d t ,d v ) i-» (d a ,du): 

= ir 2T ' ~ ^ T ' iT ' 2 + 3Ap + 1 ( 2 ^ 2 ' + 2 ^ T ^)) 

b -^T 1 = 2^(-40-# 2 T 1 '-3A p + i (2^ + 20 2 T 1 '^)). 

From the expression we notice than the coefficient of dt is of order T'^ 1 / 2 ' and 
the coefficient of d v is of order T' - ^ 1 / 2 ). This reflects the fact that the symplectic 
transformation stretches along the v direction and contracts along the t direction 
in order to put the differential in normal form. This in turn implies that the shape 
of the invariant curves is elongated in the t direction (as we can notice in figure 7). 
We are left with computing the order of magnitude of derivatives of a\ and a 2 . It 
is convenient to notice that d(F 2 ) is close to a square of a matrix, i.e. if we write 
T{ = T'(l - 5) and = T'(l + 8) we obtain: 

at = 1 + T'4> + i(4> - T'/2) a 2 = -T r <j) + + T'/2). 

ai = (aiai + a 2 a 2 ) - 20 2 T' 2 <5 2 - 2i<jy 2 T' 5 + i'<pT l2 6 2 

a 2 = (a 2 a 1 +a 2 a 1 )-24>T'5 + 24>T ,2 5 2 ~2i4> 2 T'5-i4>T' 2 S 2 

As 5 is of order 1 (and limiting to as T" — > oo). the error term is of order at most 
T', whereas the main term is of order T' 2 . Now we differentiate a\ and a 2 with 
respect to t and v: 

d t ai = <j) (T' + i) d t a 2 = '4>{-T' + i) 

d v a x = T'd t ai + T"(0 - i/2) d v a 2 = T'd t a 2 + T"(-<j> + i/2). 

In order to obtain an upper bound on such derivatives, we will consider <p n ' ~ 1 
regardless of the fact that even derivatives will be of order T' _1 . To this extent, we 
observe that all terms containing second (and higher) derivatives of T will appear 
in terms of lower order than the dominant T' for dt and T' 2 for d v . By direct 
inspection, the same statement is true for terms containing 8 in the expression for 
a, (in fact S ~ T"/T'). Therefore if we restrict to the maximum order: 

dtd l v a 1 \ max = T' l <t > {k+l+2) (T' + l ) d^d l v a 2 \ max = r l ^ k+l+2 H-T' + i). 

Now recall that we were to compute derivatives with respect to the (u, u) variables 
and as such we should recall that the coefficients in front of dt and d v are of order 
respectively T'^ 1 / 2 ^ and T'( -1 / 2 ). This means that we obtain: 

fl£4«i < r'( fc +w a+1 d h u d l u a 2 < T'( fc+/ )/ 2+1 , 

which in turn implies: 

flfi^oi % T'^' 2+2 d h u d l u a 2 < T'^' 2+2 . 
Therefore we obtain the required estimate for A^ n \ i.e.: 



which concludes the proof. 



□ 
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Figure 7: Elliptic island of period 2 and type (+); each smaller picture is an enlarged 
portion of the previous one. In the big picture we can notice the 2-periodic islands (bottom 
center and upper left) at two suitable intersection points of i and E 1 . 



